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We present a concise review of the theoretical status of the rare semileptonic B — > X a l + l~ 
decays in the standard model. Particular attention is thereby devoted to the recent theoretical 
progress concerning, on the one hand the next-to-next-to-leading order QCD calculation and, 
on the other hand the analysis of phenomenological important subleading electroweak effects. 



1 Introduction 

The rare semileptonic b — ► s£ + £~ transitions have b een ob served for the first time by Belle and 
BaBar in form of the e xclusi ve B — > K^£ + £~ modes E^lsil Recently also inclusive measurements 
have become available! 1 11 Like all other flavor-changing-neutral-current processes, these chan- 
nels are important probes of short-distance physics. Their study can yield useful complementary 
information, when confronted with the less rare b — > S7 decays, in testing the flavor sector of 
the Standard Model (SM). In particular, a precise measurement of the inclusive B — > X s £ + £~ 
mode would be welcome in view of New Physics (NP) searches, because it is amenable to a clean 
theoretical description for dilepton invariant masses, m« = q 2 , below and above the cc reso- 
nances, namely in the ranges 1 GeV 2 < q 2 < 6GeV 2 and q 2 > 14.4 GeV 2 P Each window has its 
own assets and drawbacks, related on the one hand to experimental issues, such as event rates, 
identification and detection efficiencies, on the other hand to theoretical questions, concerning 
the significance of parametric errors, perturbative and non-perturbative effects. In the following 
we present results that will cover both regions, performing a thorough study of the associated 
uncertainties. Unfortunately, however, it is not possible to compare these predictions with ex- 
periment, as it is still to early to fit a q 2 distribution to the data. To allow the comparison with 
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the existing inclusive measurements, we thus reexamine the SM prediction for the B — ► X s £ + £~ 
branching ratio, BR«, including all known QCD and electroweak effects. Finally, we update the 
SM result for the position of the zero of the Forward-Backward (FB) asymmetry, q$. 

2 Theoretical Framework 

The calculation of the partonic decay rate of B — > X s £ + £~ consists of several parts that are 
worth recalling. Perturbative QCD effects play an important role, due to the presence of large 
logarithms of the form L = lnmb/M w , that can be resummed using the machinery of Operator 
Product Expansion (OPE) and Renormalization Group (RG) improved perturbation theory. 
Factoring out the Fermi constant G^ and the electromagnetic coupling a, the amplitude receives 
contributions of 0{a™L n+1 ) at the Leading Order (LO), of 0{a™L n ) at the Next-to-Leading 
Order (NLO), and of 0{a n s L n - 1 ) at the Next-to-Next-to-Leading Order (NNLO) in QCD. 

To achieve the necessary resummation, one works in the framework of an effective low-energy 
theory with five active quarks, three active leptons, photons and gluons, obtained by integrating 
out heavy degrees of freedom characterized by a mass scale M > M w . At LO in the OPE the 
effective on-shell Lagrangian relevant for the b — > s£ + £~ transition at a scale fi is given by 

AG 

£ c s = -Cqcdxqed + -^=-V t * s V t b Ci(fi) Qi . (1) 

Here the first term is the conventional QCD-QED Lagrangian for the light SM particles. In the 
second term Vij denotes the elements of the Cabibbo-Kobayashi-Maskawa matrix and Cj(/i) are 
the Wilson coefficients of the corresponding operators Qi built out of the light fields. 

Neglecting for the moment subleading electroweak effects, as well as the QCD penguin oper- 
ators Qs~Qq, that are suppressed by small Wilson coefficients, the remaining physical operators 
arising in the SM can be written as 

Qi = (s L ^T a c L )(c L ^T a b L ) , Q 2 = (s L ^c L )(c L ^b L ) , 

Q? = jim h {s L a^b R )F^ , QI = \m b {s L a^T a b R )Gl v , (2) 

Q 9 = f^7/A)£|(W) , Qio = |r(sL7/A)£#7 M 75^) , 

where the sum over £ extends over all lepton fields, e (g) is the electromagnetic (strong) coupling 
constant, and qn are the chiral quark fields, F^ u (G^ u ) is the electromagnetic (gluonic) field 
strength tensor, and T a are the color matrices, normalized so that Tr(T a T 6 ) = 5 ab /2. 

Unlike the case of b — > 57, the b — > s£ + £~ amplitude involves large logarithms even in the 
absence of QCD interactions, since the current-current operators Q\ and Q2, as well as Q3-Q6 
mix into the vector-like semileptonic operator Qg at the one-loop level. QCD contributions 
from the magnetic operators Qj and Q|, and the axial- vector- like semileptonic operator Qiq 
enter formally at the NLO, but turn out to be numerically an 0(1) correction to the LO result. 
In order to gain an accuracy below the 10% level on the b — > s£ + £~ decay rate a complete 
resummation of NNLO QCD logarithms has thus to be p erformed. Contrary thereto, a similar 
precision is achieved for b — > S7 already at the NLO ^^accentuating once again the difference 
between rare and radiative modes in view of the RG improved perturbation theory. 

3 Recent Perturbative Standard Model Calculations 

The aforementioned NNLO QCD computation has required the computation of i) the 0(a s ) 
corrections to the relevant Wilson coefficients^ ii) the 0(a s ) contributions to the associated 



matrix elements j ^ ^ ^ ^ ^ ^ and in) the 0(a 2 ) Anomalous Dime nsion Matrix (ADM) de- 
scribing the mixing of physical dimension-five and six operators! ^ | 15 | Nearly all the ingredients 
of the NNLO QCD calculation involve a considerable degree of technical sophistication and 
have been performed independently by at least two groups, sometimes using different methods. 
However, the most complex part of the whole enterprise, the calculation of the three-loop 0(a 2 ) 
ADM describing the mixing of Qi—Qe into Q1-Q9 has been completed only very recently!^ Al- 
though the 0(a 2 ) matrix elements of Qs-Q§ have not been calculated so far, 6 the NNLO QCD 
computation of b — ► s£ + £~ can be called practically complete, as the 0(a 2 ) matrix elements of 
Qi and Q2, as well as Qgp^ evaluated for arbitrary q 2 , are now also available. In their sum 
the numerical impact of the NNLO QCD corrections on the differential decay rate amounts to 
around -20% (-25%) in the low-q 2 (high-9 2 ) region, and leads to a reduction of the renormal- 
ization scale uncertainties from around ±20% (±15%) to ±5% (±3%). In the case of the FB 
asymmetry, NNLO QCD corrections are not less important, as they shift by around ±15% 
and reduce the renormalization scale dependences from around ±20% to ±3%. 

In contrast to the track record of QCD corrections, the possible importance of subleading 
electroweak effects in b — > s£ + £~ has been realized only quite recently. As shown in the case 
of radiative decaysP^ they may be as important as the higher order QCD effects. BR^ is 
generally parameterized in terms of a, but the scale at which a should be evaluated is, in 
principle, undetermined until higher order electroweak effects are taken into account. This has 
led most authors to use ce(m«) « a(mb) 1/133. Indeed, in the absence of an 0(a) calculation, 
there is no reason to consider a(mu) more appropriate than, say, a(M w ) « 1/128. As BR^ is 
proportional to a 2 , the ensuing uncertainty of alm ost 8% is, compared to the precision achieved 
in the QCD calculation, not at all neglig 

iblePEHl 

This uncertainty, which has affected previous 
analyses, has been recently reduced to 2%, by a calculation of the numerical dominant subleading 
electroweak effects showing that after the inclusion of the latter corrections the appropriate 
prefactor in BR^ is in fact a(m&). Due to accidental cancellation, the impact of the calculated 
electroweak effects is however rather limited: it amounts to almost —2% in BR«, whereas g 2 is 
changed by around ±2%, and leaves the scale uncertainties practically unchanged. 



4 Phenomenology 



In order to cancel the strong m^-dependence of the differential B — * X s £ + £~ decay rate, it is 
customary to normalize the latter to the experimental value of the Branching Ratio (BR) for the 
inclusive semileptonic decay BR[i? — > X c £v\. However, this normalization introduces a strong 
m c -dependence, whi ch is n ot known very accurately. The ensuing uncertainty is about 8%0An 
alternative procedur J^l-^ l consists in normalizing the B — > X s £ + £~ decay width to T[B — > X u £v], 
and then to express BR[1? — ► X u £u] in terms of BR[1? — > X c £v\ and of the ratio 
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which can be computed with better accuracy.^ In other words, defining 
write the normalized differential decay rate as 
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b In B — > X s 7 the 0(o? s ) matrix elements of Q3-Q6 reduce the branching ratio by around 1% 
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+(1 + 2s) Cl s R (s) + Cfl R (§) J +12 Re {C^%(s)CI%(sr) + \ , (5) 

and the ef fective Wilson coefficients carrying the label R, include all real and virtual 0( a s ) 
corrections p^^ whereas the last term denotes the finite 0(a s ) bremsstrahlungs corrections!*^ 
The use of the 6-quark pole mass in Eq. © and in the calculation of the semileptonic width 
leads to large perturbative QCD corrections both in the numerator and denominator of Eq. @. 
Since this is an artifact of the choice of scheme for the 6-quark mass, they clearly tend to cancel 
in the ratio. As a second improvement with respect to previous analyses, we keep terms through 
0(a 2 ) in the denominator and expand the ratio in Eq. (J1J) in powers of a s . By making explicit 
the cancellation of large contributions, the convergence and stability of the perturbative series 
improves, as we have verified explicitly. Like in previous analyses, due to the peculiarity of the 
perturbative expansion for b — » sl + £~, we retain some large scheme-independent higher order 
terms in the amplitude squared: in particular, all terms in Eq. © that are quadratic in the 
effective Wilson coefficients are expanded in a s and terms up to 0(a s ) are kept. A careful study 



showsj^that with the two aforementioned improvements the residual scale dependence is about 
50% smaller than in all preciding analyses. 

After including power corrections of 0(l/m^f^ and 0(l/m 2 )^ to the differential decay 
rate of B — > X s £ + £~, as well as 0(l/m?) corrections to T[B — > X c £u] and T[B — > X u lv\, and 
expanding Eq. © in inverse powers of raj, and m c , we obtain for the normalized differential rate 
integrated over the low-g 2 region 

nn /i r> a v2 <r „2 ^ ar^Xr^ — ( 1 K7+°- n \ +0.07i +O.O61 +0.05i +0.05i \ v 1 n -6 (c\ 

BK«(lGeV <q < b GeV J - ^1.57_ . 10 | mt _ .07^-0 .07lscale-0 .oslcc-o m\c) x 10 • W 

Including other subleading parametric uncertainties, the total error is about 10% and is dom- 
inated by the uncertainty on the i-quark mass, which will soon be reduced by a factor of two 
by CDF and DO. Moreover, the substantial error from the 6-quark mass is an artifact of the 
employed scheme, which could be reduced drastically by changing the latter. Hence it should 
not be interpreted ciS ct limitation. 

Similarly, we find for the high-g 2 window 

nn („2 v. -ia a n a \ T 2\ _ (a no+0.71 1 +0.24i +0.13| +0.13 1 +0.12i \ 1n _7 /„n 

iiK u {q > 14.4 GeV J - ^4.02_ 7 J mfc _ 23 | mt _ 13 | scalc _ 13 | cE _ 012 | c J x 10 , (7) 

which is dominated by an error of around 15% related to the kinematical cut on g 2 01n fact, 
even though Eq. © is not explicitly affected by 0(1/ mi,) corrections, the physical observable 



is sensitive to 0(l/mb) terms, that arise from the well-known relation between the mass of the 
£?-meson and the 6-quark in heavy quark effective theory. 

Finally, our prediction for the integral of the partonic decay rate over the full spectrum reads 



BR« {q 2 > Am 2 ) = ( 4 - 58 ± °- 18 scalc ± 0.66 para ) x 10~ 6 , (8) 

where the parametric error takes into account, besides the power corrections of 0(l/m^) and 
0(l/m 2 ), a guesstimate of the uncertainty related to higher cc re son ances, evaluated by means 
of experimental data on e + e~ — > X c using a dispersion relation!^ The numerical size of the 
latter correction can be easily assessed from the left plot of Figure ^ c showing the NNLO QCD 
prediction for the differential branching ratio as a function of q 2 . 

Unfortunately for all of us who hoped to discover NP in rare b — > s tra nsiti ons, the estimate 
in Eq. (jHJ) compares fairly good with the recent experimental world average 1 ^ 

BR^=(6.2l^| stat i^ syst )xlO- 6 . (9) 

Even worse, it agrees amazingly well with the preliminary result by Belief which suggests a 
somewhat lower central value for the combination of all inclusive measurements. 

Facing such bad news, lets turn our attention to the FB asymmetry, which in its so-called 
normalized form is defined as 

A FB (s) = - i94 _ g _ 1 — J dcos9 e L ^ J J n ^sgn(cos^), (10) 



dT[B X s £+£~]/ds J_ x ' dsdcos( 



where 



1 dcoB0t dsdcose £ s ^ cos0 ^ = i^s {—) (i--) 

x|-6Re(cf ff FB (S)^ FB ( S T) -3fiHe(^S B («)^5 FB (fi)*) +^ ems (^)| , (11) 

and 9e denotes the angle between the momentum of the positively charged lepton and the B- 
meson in the rest frame of the lepton pair. Th e effect ive Wilson coefficients in Eq. (|1U|) take into 
account all 0(a s ) real and virtual corrections !^^' " as indic ated by the subscript FB, whereas 
the last term encodes the 0(a s ) bremsstrahlung corrections}^ which have not been included 
in our analysis of q^, since they turn out to be below 1%. The NNLO QCD prediction of 
the normalized FB asymmetry as a function of q 2 is presented in the right plot of Figure ^ 
which also illustrates the numerical size of the long-distance contributions due to intermediate 
cc resonances. 

Since q 2 is known to be especially sensitive to physics beyond the SM, a comprehens ive 
study of the residual theoretical error attached to it is important. As pointed out earlierp^ 
the renormalization scale dependence of q 2 . computed from Eq. (|llj) is rather small. However, 
an alternative way to estimate the residual theoretical error follows from the observation that 
should be independent of the normalization for the FB asymmetry. As both numerator and 
denominator of Eq. (fTU|) are truncated series in a s , one can expand Eq. (fTU|) in a s , making q$ 
sensitive to different combinations of higher order terms, that depend on the adopted normal- 
ization. A thorough analysis showsp^ that the dependence of q$ on the specific method used 
to compute it is generally larger than the scale uncertainty associated with it. Implementing all 
known perturbative and non-perturbative corrections, and assigning to the central value a rather 
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conservative theoretical error of 6%, that covers the whole range of possible values obtained by 
changing the normalization in Eq. (jlOj) . we finally get 

ql = (3.76 ± 0.22 theory ± 0.24 m6 ) GeV 2 , (12) 

with a total error close to 9%. Also in this case the error due to the 6-quark mass can be 
drastically reduced by performing the calculation in a different mass scheme. 



5 Outlook 



Rare semileptonic B decays are going to play an important role in the search for NP at the 
£?-factories and upcoming flavor physics experiments at the Tevatron and the LHC. From the 
theoretical point of view the inclusive mode in both the low-t/ 2 and high-g 2 window is particular 
interesting, since it can be accurately computed in the SM. Given this situation, a measurement 
of the differential decay rate integrated separately over the two reference regions would be 
desirable. Furthermore, as the FB asymmetry can change drastically in scenarios of NP, even a 
rather crude measurement of its shape would either rule out large parts of the parameter space 
of the underlying model or show clear evidence for physics beyond the SM. 
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